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INTRODUCTION 


The description of the computer program is designed to display and evaluate the 
comparison of multilayer, dielectric, thin-film optical filter transmittance values with 
corresponding wavelength data for normal incidence only. 

The program is written for the PDP-5 and PDP-8 combined computer system 
which has an intercommunication buffer, a 30 -G cathode ray tube (CRT) display, a deck 
tape system, and teletype and punched-tape input and output equipment. 

This technical memorandum is a supplement to NASA Technical Memoran- 
dum X-58013, November 1967. 


COMPUTER PROGRAM 


The "Thin-Film Filter Computer Program System" has been written to operate 
with the PDP-5/8 computers linked with an interface connection. The system uses the 
following peripheral equipment: a teletype input -output system and a 30-G display unit 
with 14 -inch CRT. 

The program system is in three parts. The main program is a FORTRAN pro- 
gram that reads the input data tape, performs the necessary calculations, sends a 
signal to the CRT link routine, and prints complete results on the teletype unit. The 
CRT display program is an assembler language program that displays the FORTRAN 
program results on the CRT; the link program is an assembler language program that 
transmits the FORTRAN program output to the PDP-5. 

The input data for a run are arranged on paper tape according to the format in 
table I. Definitions of the variables listed in table I are found in table II. 


♦Programing by William J. Fehrenkamp and Marian VanZant; program analysis 
by James W. Martin. 




The data for adjusting the thicknesses may be prepared on paper tape or typed 
manually on the teletype as the program proceeds. A sample data tape is shown in 
table III. The data on the paper tape are transmitted to the computer by the teletype 
unit. 


Output on the teletype unit is shown in tables IV to VII. An output sample for CRT 
can be found on pages 6 and 8 of NASA TM X-58013 entitled "Evaluation of Optical Thin- 
Film Filters." Following table VII are the listings and flow charts (figs. 1 and 2) for 
the three computer programs. 

The loading sequence is in two main parts. The first part is for the PDP-5, and 
the second part is for the PDP-8. 

1. PDP-5 — Load CRT display routine and floating point package A in either 
order with the binary loader; branch to 0400 to start the display routine. 

2. PDP-8 — Load according to the following order: 

a. FORTRAN operating system no. 1 with binary loader 

b. Assembler language link routine with binary loader 

c. FORTRAN program object tape with low -speed reader option 
Branch to 0200 and begin. 
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TABLE I. - INPUT SEQUENCE AND FORMAT 


Input sequence 

Format 

N 

I 

M 

I 

XNO 

E 

XNS 

E 

XN X 

E 

XN 

E 

n 


XLAM 2 

E 

XLAM 

E 

m 


THj 

E 

TH 

E 

n 


1 

NINP, I, TH. 

I, I, E 

a 

1 1 

I, I, E 


' Repeat until NINP = 0 



a 


Repeat for each altered set of thicknesses. 







TABLE II. - INPUT DEFINITIONS 


Initial 

Definition 

N 

Number of layers 

M 

Number Of lambdas 

XNO 

Incident medium refractive index (N q ) 

XNS 

Substrate refractive index ^N g ) 

XNi 

Index of refraction of layer 

XLAMi 

Wavelength of light (lambda) 

THi 

Thickness of a layer 

Thickness Alterations 9, 

NINP 

If NINP - 0, last thickness change 
If NINP = 1, more changes coming 

I 

That layer where thickness is to be changed 

Th(i) 

The changed thickness 


£1 ‘ 

This section is repeated every time a new altered thick- 
ness table is desired. 
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TABLE III. - INPUT TAPE FOR INITIAL SET OF THICKNESSES 2- 


Input parameters 

Wavelength A 

Thickness t 

18 

. 40E-04 

.643E-05 

15 

. 41E-04 

.751E-05 

1.0 

. 42E-04 

.890E-05 

1.52 

. 43E-04 

•751E-05 

1. 38 

.44E-04 

.890E-05 

2.30 

. 45E-04 

•751E-05 

1.38 

. 46E-04 

.890E-05 

2.30 

. 47E-04 

•751E-05 

1.38 

. 48E-04 

.890E-05 

2.30 

. 49E-04 

•751E-05 

1.38 

, 50E-04 

.890E-05 

2.30 

.51E-04 

•751E-05 

1.38 

. 52E-04 

.890E-05 

2.30 

. 53E-04 

•751E-05 

1.38 

. 54E-04 

.890E-05 

2. 30 


•751E-05 

1.38 


.890E-05 

2.30 


•751E-05 

1.38 



2.30 



1.38 



2.30 




^nput order: column 1 followed by column 2 followed by column 3. 
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TABLE IV. - OUTPUT OF INITIAL SET PARAMETERS 


NUMBER OF LAYERS = l8 
NUMBER OF LAMBDAS = 15 

XNO = 1.0 


XNS = 1.52 

LAYER INDEX OF REFRACTION 

+1 1.38 

+2 2.30 

+3 1.38 

+4 2.30 

+5 1.38 

+6 2.30 

+7 1-38 

+8 2.30 

+9 1.38 

+10 2.30 

+11 1.38 
+12 2.30 

+13 1.38 

+14 2.30 

+15 1. 38 

+l6 2.30 

+17 1.38 

+18 2.30 

LAMBDA 

. 40E-04 
. 41E-04 
. 42E-04 
.43E-04 
. 44E-04 
.45E-04 
.46E-04 
. 4-7E-04 
.48E-04 
.49E-04 
. 50E-04 
. 51E-04 
. 52E-04 
. 53E-04 
.54E-04 
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TABLE V. - OUTPUT OF INITIAL SET, \ VERSUS T 


LAYER 

THICKNESS 

OF LAYER 

+1 

.643E-05 


+2 

.751E-05 


+3 

.890E-05 


+4 

.751E-05 


+5 

. 890E-05 


+6 

.751E-05 


+7 

.890E-05 


+8 

.751E-05 


+9 

.890E-05 


+10 

•751E-05 


+11 

. 890E-05 


+12 

.751E-05 


+13 

.890E-05 


+l4 

.751E-05 


+15 

.890E-05 


+l 6 

. 751E-05 


+17 

.890E-05 


+18 

. 751E-05 


LAMBDA 


TRANSMITTANCE 

+0 . 399999E-4 

+0.987624E+0 

+0.409999E-4 

+0.942006E+0 

+0 . 419999E-4 

+0.972206E+0 

+0 . 429999E-4 

+0.992719E+0 

+0.439999E-4 

+0.910667E+0 

+0.1+1+9999E-4 

+0.892187E+0 

+ 0 . 1 + 59999 E— L 

+0.984376E+0 

+0.469999E-4 

+0 . 898954E+0 

+0.1+80000E-1+ 

+0.693483E+0 

+0 ... 489999E-4 

+0.785219E+0 

+0.1+99999E-4 

+0.779257E+0 

+0.509999E-4 

+0.117121E+0 

+0.519999E-4 

+0 . 248735E-1 

+0.529999E-4 

+0.829176E-2 

+0.539999E-4 

+0.369789E-2 
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TABLE VI. - OUTPUT OF ADJUSTED SET 


1 12 659E-05 

Thickness adjustments 

0 13 76TE-05 


LAYER THICKNESS 

+1 +0.642998E-5 
+2 +0.750999E-5 
+3 +0.889999E-5 
+k +0.750999E-5 
+5 +0.889999E-5 
+6 +0.750999E-5 
+7 +0.889999E-5 
+8 +0.750999E-5 
+9 +0.889999E-5 
+10 +0.750999E-5 
+11 +0.889999E-5 
+12 +0.658999E-2 
+13 +0.767000E-2 
+l4 +0.750999E-5 
+15 +0.889999E-5 
+16 +0.750999E-5 
+17 +0.889999E-5 
+18 +0.750999E-5 


LAMBDA 

+0.399999E-U 
+0-L09 999E-4 
+0 . 419999E-4 
+0.1+29999E-4 
+0. 439999E-1+ 
+0.4U9999E-4 
+0.459999E-4 
+0. L69999E-J+ 
+0.480000E-4 
+0.489999E-4 
+0 . 1 +99999E— 1+ 
+0.509999E-4 
+0.519999E-4 
+0.529999E-U 
+0.539999E-4 


TRANSMITTANCE 

+0.832921E+0 
+0.789471E+0 
+0.984829E+0 
+0.935955E+0 
+0.937062E+0 
+0.887739E+0 
+0 . 870U50E+0 
+0 . 682983E+0 
+0.389687E+0 
+0 . 206427E+0 
+0 . 374616E-1 
+0.278002E-1 
+0.591010E-1 
+0.127748E+0 
+0.388402E-1 



TABLE VII. - OUTPUT OF ADJUSTED SET 


01 8 .699E-05 

01 9 •777E-05 Thickness adjustments 2 

00 10 .688E-05 

LAYER THICKNESS 

+1 +0.61*2998E-5 

+2 +0.750999E-5 

+3 +0.889999E-5 

+4 +0.750999E-5 

+5 +0.889999E-5 

+6 +0.750999E-5 

+7 +0.889999E-5 

+8 +0.698999E-5 

+9 +0 776999E-5 

+10 +0.687999E-5 

+11 +0.889999E-5 

+12 +0.658999E-2 

+13 +0.767000E-2 

+14 +0.750999E-5 

+15 +0.889999E-5 

+16 +0.750999E-5 

+17 +0.889999E-5 

+18 +0.750999E-5 

LAMBDA TRANSMITTANCE 

+0.399999E-4 +0.944412E+0 

+0 . 409999E-4 +0 . 620119E+0 

+0 . 419999E-4 +0.869817E+0 

+0.429999E-4 +0.796673E+0 

+0. 439999E-4 +0.866307E+0 

+0.449999E-4 +0 . 848100E+0 

+0.4-59999E-4 +0.71T866E+0 

+0 . 469999E-4 +0 . 386330E+0 

+0.480000E-4 +0.103591E+0 

+0.489999E-4 +0.481866E-1 

+0.499999E-4 +0.153897E-1 

+0 • 509999E-4 +0 . 173737E-1 

+0 • 519999E-4 +0 . 177824E-1 

+0.529999E-4 +0.583763E+0 

+0.539999E-4 +0.664485E-1 
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o o o 


MAIN PROGRAM (FORTRAN) 


C; OPTICAL THIN-FILM FILTER EVALUATION 
C; MANNED SPACECRAFT CENTER 
C; HOUSTON TEXAS 
C; 

C; READ INITIAL PARAMETERS ON TELETYPE 
C; (MANUAL OR PAPER TAPE) 

C; 

DIMENSION XN(25) ,TH(25) ,XLAM(25) ,T(25) 

100; FORMAT ( I ) 

200; FORMAT (E) 

TYPE 101 

101; FORMAT ( / ," NUMBER OF LAYERS = ") 

ACCEPT 100 ,N 
TYPE 102 

102; FORMAT ( / /'NUMBER OF LAMBDAS = ") 

ACCEPT 100 ,M 
TYPE 103 

103; FORMAT (/,"XNO = ") 

ACCEPT 200, XNO 
TYPE 104 

104; FORMAT (/,"XNS = ") 

ACCEPT 200 ,XNS 
TYPE 105 

105; FORMAT ( / , "LAYER ", "INDEX OF REFRACTION”,/) 
DO 20 1=1 ,N 
TYPE 106,1 
106; FORMAT (i) 

ACCEPT 200 ,XN(l) 

20; CONTINUE 
TYPE 107 

107; FORMAT(/, "LAMBDA" ,/) 

DO 30 1=1 ,M 
ACCEPT 200,XLAM(I) 

30; CONTINUE 

; INPUT FIRST SET OF THICKNESSES 
TYPE 108 

108; FORMAT (//'LAYER "/'THICKNESS OF LAYER",/) 
DO 40 1=1 ,N 
TYPE 109,1 
109; FORMAT (I) 

ACCEPT 200,TH(I) 

40; CONTINUE 
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MAIN PROGRAM (FORTRAN) - Continued 


; T » s COMPUTED IN THIS SECTION 

A 


COMPUTE DEL AND M MATRIX 
998; DO 999 ILAM=1,M 

DEL=( ( 6 . 2831853 ) *XN ( 1 ) *TH( 1 ) ) /XLAM(IIiAM) 

A=COSF(DEL) 

B=SINF(DEL)/XN(l) 

C=XN ( 1 ) *SINF ( DEL ) 

D=A 

IF(N-l)2,4,2 

; HERE FOR N>1 COMPUTE M =M 

n 

2; DO 3 1=2 ,N 

DEL“{ ( 6 . 2831853 ) *XN ( I ) *TH( I ) ) /XLAM( ILAM) 

ai=cosf(del) 

bi=sinf(del)/xn(i) 

■ C1=XN(I)*SINF(DEL) 

D1=A1 

A2=A*A1-B*C1 

B2=A*B1+B*D1 

C2=C*A1+D*C1 

D2=D*D1-C*B1 

A=A2 

B=B2 

C=C2 

D=D2 

3; CONTINUE 
; T EQUATION 

A 

4 ;T( ILAM) =4 .0/(2. 0+A*A*XNO/XNS+D*D*XNS/XNO+C*C/ ( XNO*XNS ) +B*B*XNO*XNS ) 
999; CONTINUE 


12 



o o o a o a a o o 


MAIN PROGRAM (FORTRAN) - Concluded 


; GO TO DISPLAY ROUTINE AND RETURN 
TYPE 110 

110; FORMAT ( / , "LAMBDA "."TRANSMITTANCE",/) 

PAUSE 3^56 

• PRINT X x , T x TABLE 

DO 60 1=1 ,M 

TYPE 300,XLAM(I),T(I) 

300; FORMAT (/,E,E) 

60; CONTINUE 

; INPUT ALTERATIONS IN THICKNESS TABLE 

70; ACCEPT 400 jNINP , I , TH ( I ) 

400; FORMAT (I, I, E) 

C; IF NOS OF INPUTS = 0, START COMPUTING 

C; IF NOS OF INPUTS IS NON-ZERO, CONTINUE TO INPUT 

IF (NINP)70,80,70 

80; TYPE 111 

111; FORMAT ( / , "LAYER "."THICKNESS",/) 

DO 90 1=1 ,N 
TYPE 500,I,TH(I) 

500; FORMAT(/,I,E) 

90; CONTINUE 
GO TO 998 
END 
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LINK PROGRAM (ASSEMBLER LANGUAGE) 


/THIS SUBROUTINE TRANSFERS BILBAO PROGRAM DATA 

/FROM THE PDP-8 TO THE PDP-5 

*6600 

SUBR, 0 



CLA 

TAD BADR 

/SET UP ADDRESS POINTER 


DCA GADR 

/ 


tad i nlam 

/GET NUMBER OF LAMBDAS AND 


DCA TEMPI 

/TRANSMITTANCES 


TAD TEMPI 

/MULT BY -3 


TAD TEMPI 

/..... 


TAD TEMPI 

/,.,.. 


CIA 

/.,... 


DCA TEMP2 

/PUT IN TEMP2 TO USE AS A COUNTER 


TAD TEMPI 

/GET NUM, OF L&T 


6716 

/TRANSFER NOM. TO. 5 

RXWORD, 

CLA 



TAD I CADR 

/GET TRANSMITTANCE DATA WORD 


6711 

/WAIT UNTIL 5 IS READY 


SKP 

/ 


JMP .-2 

/ 


6716 

/TRANSFER TO 5 


ISZ CADR 

/INCREMENT POINTER 


ISZ TEMP2 

/INCREMENT COUNTER AND SKIP IF FINISHED 


JMP EXWORD 

/NOT FINISHED, GET ANOTHER DATA WORD 


JMP I SUBR 

/FINISHED, EXIT 

TEMPI , 

0 


TEMP2, 

0 


NLAM, 

704? 


BADR, 

7051 


CADR, 

7051 




i)72r\ 

i«5 u 

READY. 

V? X 


BEGIN 
SA 0400 


GO TO 
GETDAT 


CLEAR AC 
LOAD AC 0477 
STORE IN 17 P.0, 


LOAD AC 
SW REGISTER 


BRIGHTEN CRT 
CLR X BFR 


LOAD AC 
WITH -20. 
STR IN CNT 


CLEAR Y BFR 
LOAD AC XPNT 
XEC X COORD 


ADD 32.T0 AC 
STORE IN XPNT 
LOAD AC (IND) 
17 P.0. 


CNT>0? 


XEC Y COORD 


XEC Y COORD 


CLEAR AC 
INCRMENT CNT 


LOAD AC 
WITH 50. 
INCREMENT 
CNT 


CNT>0? 


Figure 2. - Cathode ray tube flow diagram. 
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NO 


2. - Continued. 








^ PROC ^ 




CLEAR AC AND LINK 


1 




GO TO (IND) 

0007 P.0. 

L > 




LOAD FLT AC [lOOO] 
NORMALIZE 
STORE IN 1000 


CLEAR AC 

LOAD AC 600 OR 384. 
STORE IN ADDR 


1 LOAD AC 477 

i 

OR 319. 
IN 17 P 

— i 

STORE 

.0. 

1 

1 

1 

i 

GO TO 
ERR2 


LOAD AC NUMBY 
STORE IN CNT 


GO TO (IND) 
0007 P.0. 






LOAD FLT 
MULT. BY 

AC ADDR 

1000] 

✓ 

t i 




Figure 2. - Continued. 
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CLEAR LINK 
LOAD AC (45) 








CRT DISPLAY ROUTINE FOR PDP-5 (ASSEMBLER LANGUAGE) 




*7 

0007 

5600 

*63 

006 3 

0000 

FIX, 

0064 

7200 


0065 

1044 


0066 

7540 


006 7 

5072 


0070 

7200 


0071 

5111 


0072 

1112 


0073 

7450 


0074 

5110 


0075 

7500 


0076 

5513 


0077 

3044 


0100 

7100 

AGN, 

0101 

1045 


0102 

7510 


0103 

7020 


0104 

7010 


0105 

3045 


0106 

2044 


0107 

5100 


0110 

1045 

FIN, 

0111 

5463 


0112 

7765 

M13, 

0113 

1000 

ERRl, 

0114 

1001 

ERR2, 


5600 

0 

CLA 

TAD 44 
SZA SMA 
JMP .+3 
CLA 

JMP FIN +1 
TAD M13 
SNA 

JMP FIN 
SMA 

JMP I ERR1 
DCA J+L 
CLL 

TAD 45 
SPA 
CML 
BAR 
DCA 45 

isz 44 

JMP AGN 
TAD 45 
JMP I FIX 

0-13 

1000 

1001 
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CRT DISPLAY ROUTINE FOR PDP-5 (ASSEMBLER LANGUAGE) - Continued 


0-400 

5201 

*400 

BEGIN, 

JMP .+1 

0401 

6721 


6721 

0402 

5246 


JMP GETDAT 

0403 

7604 


LAS 

0404 

6074 


607^ 

0405 

6051 


6051 

0406 

7200 


CLA 

0407 

1237 


TAD C20 

0410 

3236 


DCA CNT 

04ll 

60 67 


60 67 

0412 

1240 


TAD C50 

0413 

2236 


ISZ CNT 

04l4 

5211 


JMP .-3 

0415 

7200 


CLA 

04l6 

1241 


TAD C477 

0417 

3017 


DCA Z 17 

0420 

1245 


TAD NUMBY 

0421 

3236 


DCA CNT 

0422 

3244 


DCA XPNT 

0423 

6061 

NXTX, 

6061 

0424 

1244 


TAD XPNT 

0425 

6057 


6057 

0426 

1242 


TAD C40 

0427 

3244 


DCA XPNT 

0430 

1417 


TAD I Z 17 

0431 

6067 


6067 

0432 

7200 


CLA 

0433 

2236 


ISZ CNT 

0434 

5223 


JMP ..-11 

0435 

5201 


JMP BEGIN +1 

0436 

0000 

CNT, 

C20, 

0 

0437 

7754 

775^ 

0440 

006 2 

C20, 

0062 

0441 

0477 

C477, 

0^77 

0442 

0040 

C40, 

40 

0443 

0577 

C577, 

577 

0444 

0000 

XPNT, 

0 

0445 

7740 

NUMBY, 

77^0 

0446 

7300 

GETDAT, 

CLA CLL 

0447 

7020 


CML 

0450 

1243 


TAD C577 

0451 

3016 


DCA Z 16 

0452 

6722 


6722 

0453 

7041 


CIA 

0454 

3245 


DCA NUMBY 
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CRT DISPLAY PROGRAM FOR PDP-5 (ASSEMBLER LANGUAGE) - Continued 


0455 

6724 

NXWORD, 

6724 

01+56 

7000 


NOP 

0457 

7000 


NOP 

0460 

7000 


NOP 

046l 

7000 


NOP 

0462 

7000 


NOP 

0463 

7000 


NOP 

0464 

6721 


6721 

0465 

7410 


SKP 

0466 

5342 


JMP PROC 

0467 

6722 


6722 

0470 

3416 


DCA I Z 16 

0471 

5255 


JMP NXWORD 

0472 

0012 

C1000 , 

0012 

0473 

3720 


3720 

0474 

0000 


0000 

0475 

0600 

ADDR, 

600 

0476 

0600 

C600, 

600 



*500 


0500 

0000 


0 

0501 

0000 


0 

0502 

0000 


0 

0503 

0000 


0 

0504 

0000 


0 

0505 

0000 


0 

0506 

0000 


0 

0507 

0000 


0 

0510 

0000 


0 

0511 

0000 


0 

0512 

0000 


0 

0513 

0000 


0 

0514 

0000 


0 

0515 

0000 


0 

0516 

0000 


0 

0517 

0000 


0 

0520 

0000 


0 

0521 

0000 


0 

0522 

0000 


0 

0523 

0000 


0 

0524 

0000 


0 

0525 

0000 


0 

0526 

0000 


0 

0527 

0000 


0 

0530 

0000 


0 

0531 

0000 


0 

0532 

0000 


0 

0533 

0000 


0 

0534 

0000 


0 

0535 

0000 


0 

0536 

0000 


0 

0537 

0000 


0 

0540 

0000 


0 

0541 

0000 


0 


21 



CRT DISPLAY PROGRAM FOR PDP-5 (ASSEMBLER LANGUAGE) - Concluded 


*5*2 


05*2 

7300 PROC, 

CLA CLL 

05*3 

4407 

JMS I Z 7 

054* 

5272 

FGET C1000 

0545 

7000 

FNOR 

0546 

6272 

FPUT C1000 

05*7 

0000 

FEXT 

0550 

7200 

CLA 

0551 

1276 

TAD C600 

0552 

3275 

DCA ADDR 

0553 

1241 

TAD 0^77 

0554 

3017 

DCA Z 17 

0555 

1245 

TAD MJMBY 

0556 

3236 

DCA CNT 

0557 

4407 

JMS I Z 7 

0560 

5675 

FGET I ADDR 

0561 

3272 

FMPY Cl 000 

0562 

0000 

FEXT 

0563 

4063 

JMS Z FIX 

0564 

7510 

SPA 

0565 

5514 

JMP. I Z ERR2 

0566 

3417 

DCA I Z 17 

0567 

2275 

ISZ ADDR 

0570 

2275 

ISZ ADDR 

0571 

2275 

ISZ ADDR 

0572 

2236 

ISZ CNT 

0573 

5357 

JMP ..-ill 

0574 

5203 

*1000 

JMP BEGIN +3 

1000 

7402 

HLT 

1001 

7402 

HLT 

1002 

5201 

JMP .-1 

ADDR 

0475 


AGN 

0100 


BEGIN 

0400 


GNT 

0436 


€1000 

0472 


C20 

0437 


C40 

0442 


0477 

0441 


050 

0440 


C577 

0443 


C600 

0476 


ERR1 

0113 


ERR2 

0114 


FIN 

0110 


FIX 

006 3 


GETDAT 0446 


M13 

0112 


NUMBY 

0445 


IKTX 

0423 


NXWORD 0455 


PROC 

0542 


XPNT 

0444 
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